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THE  EFFECTS  OF  FIELD  ERRORS  ON 
LOW-GAIN  FREE  ELECTRON  LASERS 

I.  Introduction 

Intrinsic  magnetic  field  errors  SB  are  present  in  any  realistic  wiggler  magnet.  Such 
errors  are  unavoidable  and  arise  from  imperfections  in  the  fabrication  and  assembly  of 
wiggler  magnets.  State-of-the-art  wiggler  construction  techniques  yield  rms  field  errors 
on  the  order^  {^B / —  0-1  —  0.5%.  These  field  errors  perturb  the  electron  beam  as  it 
propagates  through  the  wiggler^"^*  and  lead  to  i)  a  random  walk  of  the  beam  centroid,  Sx, 
ii)  variations  in  the  axial  beam  energy,  ^7*,  and  iii)  variations  in  the  relative  phase  of  the 
electrons  in  the  ponderomotive  potential,  If  left  uncorrected,  field  errors  ultimately 
decrease  free  electron  laser  (FEL)  gain^“^®  (this  reduction  becomes  more  significant  for 
long  wigglers).  Reduction  in  gain  may  occur  from  a  loss  of  overlap  between  the  radiation 
and  electron  beam  (due  to  large  Sx)  or  from  a  loss  of  FEL  resonance  (due  to  large  Sr/>). 

The  initial  research  on  the  effects  of  field  errors,  for  the  most  part,  was  primarily 
concerned  with  the  random  walk  Sx.  It  has  been  shown  that  the  random  walk  Sx  may  be 
effectively  controlled  by  i)  transverse  beam  focusing®”®  (finite  k^,  where  kp  is  the  betatron 
wavenumber)  and  by  ii)  periodic  beam  steering.®”^®  By  using  either  one  or  a  combination 
of  beam  focusing  and  periodic  steering,  in  principle,  the  random  walk  Sx  may  be  kept  as 
small  as  desired.  The  major  conclusions  of  the  present  work  are  the  following.  Given  that 
the  random  walk  Sx  may  be  effectively  controlled,  the  phase  deviation  Stp  is  the  primary 
physical  parameter  characterizing  loss  of  gain  for  FELs  in  the  low-gain  regime.®”^®  In 
particular,  in  order  to  avoid  significant  reduction  in  gain,  it  is  necessary  that  \Srj}\  «  27r. 
In  addition,  transverse  beam  focusing  is  not  effective  in  controlling  Specifically,  it  may 
be  shown  that  the  mean  phaise  deviation  {Sif})  is  independent  of  transverse  focusing  (inde¬ 
pendent  of  kp),  where  (...)  signifies  an  ensemble  average.  Furthermore,  beam  steering®”*® 
may  be  used  to  reduce  \St}}\  when®  Ls  <  where  Ls  is  the  length  over  which  the  steer¬ 
ing  is  performed  and  X/j  =  2irfkp.  As  an  example,  for  k^  =  0  and  one  steering  segment, 
{St{})  =  (1/3)(6'0)n>  where  {SiI))n  is  the  value  in  the  absence  of  steering. 

As  a  further  motivation,  it  is  appropriate  to  consider  some  aspects  of  wiggler  design. 
Typically,  when  “ordering”  a  wiggler  from  a  vendor,  limits  are  placed  on  SBrmi  =  {SB^)^^^ 
and  I  f  dzSB\.  To  meet  these  specifications,  the  vendor  may  “arrange”  the  magnet  poles 
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(i.e.,  actual  pole  rearrangements,  the  use  of  shims,  judicious  magnet  selection,  etc.)  in 
an  optimum  sequence^^"^®  such  that  |  J  dz6B\  is  minimized.  This  and  other  research^"^* 
indicates,  however,  that  for  low-gain  FELs  the  optimum  “figure  of  merit”  to  minimize  is 
not  the  line  integral  |  f  dz6B\,  but  the  magnitude  of  the  phase  deviation  \6rl}\. 

In  the  following,  the  effects  of  random  transverse  magnetic  field  errors,  6Bx{^),  on 
the  performance  of  low-gain  FELs  are  studied  analytically  and  numerically.  In  particular, 
the  transverse  displacement,  parallel  energy  variation  and  relative  phase  deviation  of  an 
electron  beam  propagating  through  a  wiggler  are  calculated  neglecting  the  effects  of  finite 
beam  emittance,  initial  beam  energy  spread  and  wiggler  field  tapering.  Furthermore,  the 
FEL  gain  in  the  low-gain  regime  is  determined  in  the  ID  limit,  assuming  a  plane  wave, 
non-diffracting  radiation  field.  Expressions  are  derived  for  a  particular  FEL  quantity  Q  for 
a  single  wiggler  realization  (a  specific  set  of  field  errors)  and  for  an  ensemble  of  statistically 
identical  wigglers.  The  ensemble  averages  (the  mean  and  the  variance)  of  the  quantity  Q 
are  useful  for  determining  the  most  probable  range  of  Q  for  a  particular  member  of  the 
ensemble.  The  remainder  of  this  paper  is  organized  as  follows.  The  random  walk  of  the 
beam  centroid  and  the  consequent  variations  in  the  axial  beam  energy  are  discussed  in 
Sections  II  and  III,  respectively.  The  deviations  in  the  relative  phase  resulting  from  the 
field  errors  are  examined  in  Section  IV.  In  Section  V,  the  effect  of  the  field  errors  on  the 
FEL  gain  in  the  low-gain  regime  is  determined.  The  benefits  of  beam  steering  are  analyzed 
in  Section  VI  and  addition  methods  for  reducing  the  detrimental  effects  of  field  errors  are 
discussed  in  Section  VII.  This  paper  concludes  with  a  discussion  and  summary  in  Section 
VIII. 
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II.  Transverse  Beam  Centroid  Deviations 

As  the  electron  beam  propagates  through  the  wiggler,  the  electrons  experience  random 
velocity  kicks  ^wjl  via  the  v*  x  8B±  random  force.  The  transverse  centroid  motion  of  an 
electron  beam  passing  through  a  wiggler  with  transverse  gradients  (weak  focusing)  and 
finite  field  errors  is  characterized  by  an  equation  of  the  form* 

d^8x/dz^  =  -k08x  +  k^a^8By/j,  (1) 

where  k.^  is  the  wiggler  wavenumber,  k^  =  kwa^/{\/2j)  is  the  betatron  wavenumber, 
^By  =  8ByfB^  is  the  normalized  field  error,  B^  is  the  ideal  wiggler  peak  magnetic  field, 
dv,  =  eByi/k^mc^ ,  7  is  the  relativistic  factor  of  the  electron  beam  and  z  is  the  axial 
propagation  distance.  The  first  term  on  the  right  represents  the  focusing  force  due  to  the 
transverse  gradients  in  the  wiggler  field,  whereas  the  second  term  on  the  right  represents 
the  random  force  due  to  the  field  errors.  This  equation  may  be  solved  to  give  the  random 
centroid  motion* 

8fif  =  ~  —  f  dz'  cos k/3{z'  -  z)8By{z'),  (2) 

7  Jo 

8x  =  — ^  f  dz' s\Ti  k0[z’ —  2)6 By [z'), 

7*0  Jo 

where  5/3*  =  5v*/c  is  the  normalized  transverse  velocity  deviation. 

Given  the  precise  functional  dependence  of  the  wiggler  errors  8By{z)  for  a  given  wig¬ 
gler,  the  above  expressions  may  be  used  to  calculate  the  transverse  orbit  deviations  5/3* (z) 
and  8x(z)  for  that  specific  wiggler.  However,  one  does  not  always  know  ahead  of  time 
the  full  functional  dependence  of  8By{z).  Instead,  one  may  know  only  certain  statistical 
properties  of  the  field  errors,  such  as  the  rms  value  8Brm.f  Hence,  it  is  useful  to  consider 
an  ensemble  of  statistically  identical  wigglers  for  which  the  statistical  properties  of  the 
field  errors  are  known.  By  performing  appropriate  averages  over  this  ensemble,  one  may 
determine  the  mean  {Q)  and  variance  <r  for  a  quantity  Q  and,  hence,  determine  the  most 
probable  range  of  a  single  realization  of  Q.  Here  and  throughout  the  following,  8By(z)  is 
assumed*  to  be  a  random,  homogeneous  function  with  zero  mean,  finite  variance  and  with 
an  autocorrelation  distance  given  by  z^y  (z^y  —  A„/2  is  assumed).  Details  of  the  statistical 
properties  of  of  the  field  errors  8By{z)  are  discussed  in  Ref.  6. 
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Statistically  averaging  over  an  ensemble  of  wigglers,  it  is  possible  to  determine  the 
mean-square  centroid  motion*  (neglecting  the  effects  of  finite  beam  emittance) 


!•> 

where  D  =  a^klj{6By)zcyl{2'y^).  Physically,  the  centroid  orbits  6x  and  represent 
diffusing  betatron  orbits  characterized  by  a  diffusion  coefficient  D.  Notice  that  by  in¬ 
creasing  by  additional  external  focusing,  one  may,  in  principle,  keep  Sxrm*  as  small 
as  desired.  (The  minimum  centroid  displacement  is  limited  by  finite  emittance  effects, 
which  are  neglected  in  the  present  discussion.)  Furthermore,  notice  that  in  the  ID  limit, 
(2fc^2)^  <<  1,  (6^1)  =  2D2  and  (Sx^)  =  2Dz^ /3,  as  found  previously  by  Kincaid.^  Hence, 
weak  focusing  (finite  k^)  is  effective  in  reducing  the  asymptotic  scaling  of  the  random  walk 
6x 

rm*  from  2*/^  to  2^^^.  To  avoid  loss  of  the  overlap  between  the  radiation  and  electron 
beam,  it  is  desirable  to  keep  {8  x^)  «  rj,  where  r,  is  the  radiation  spot  size.  A  detailed 
discussion  of  transverse  orbit  deviations  arising  from  random  field  errors  in  various  wiggler 
configurations  is  given  in  Ref.  6  and  summarized  in  Appendix  A. 
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III.  Axial  Beam  Energy  Variations 

Not  only  do  the  field  errors  perturb  the  transverse  motion  of  the  electrons,  they  also 
perturb  the  axial  motion.  This  is  true  since  a  static  magnetic  field  conserves  total  electron 
energy.  The  axisd  motion  may  easily  be  calculated*  using  the  above  expressions  for  the 
transverse  motion  along  with  =constant.  One  may  calculate  various  statistical 

moments  of  the  axicil  motion,  such  as  the  mean  axial  energy  variation  {S’^z)  =  (7z)  —  7z0- 
For  example,  the  mean  energy  variation  for  a  helical  wiggler  with  transverse  focusing  is 
given  by 

{h^)  _  (1  +“L/4).2  ,.2  ^ 


where  the  limit  [2k0z)^  >>  1  has  been  assumed  and  Zc  =  Zci  =  Zcy  A  detailed  discussion 
of  the  axial  energy  variation  for  various  wiggler  configurations  is  given  in  Ref.  6. 


Statistically,  (S'Yz)  may  be  interpreted  as  an  effective  energy  spread  due  to  field 
errors.®  This  effective  energy  spread  may  lead  to  a  loss  of  FEL  resonance.  Heuristi- 
cally,  in  order  to  maintain  resonance,  one  expects  that  in  the  low  or  high  gain  regime 
the  effective  energy  spread  must  be  small  compared  to  the  intrinsic  FEL  efficiency  t}, 
I  {^7*)  1/7*0  <  In  the  trapped  particle  regime,  maintaining  resonance  implies  that  the 
effective  energy  spread  must  be  small  compared  to  the  depth  of  the  ponderomotive  well, 
1(^7*)  1/7*0  <  |e$p|/(7Tnc*),  where  is  the  ponderomotive  potential.  For  example,  in  the 
low-gain  regime,  t]  =  \/{2N),  where  N  is  the  number  of  wiggler  periods.  The  inequality 
l(^7*)l/7*o  <  V  implies  6Brm»  <  l/C’riV)  ~  0.3%  for  N  =  100  (where  ~  1  has  been 
assumed  and  8Brm»  =  {8Bx)^^^  =  {SBy)^^^).  One  should  keep  in  mind  that  other  sources 
exist  which  contribute  to  the  effective  energy  spread  (e.g.,  initial  beam  emittance  and  en¬ 
ergy  spreads  due  to  transverse  gradients  in  the  wiggler  fields)  and  these  factors  should  be 
considered  in  a  complete  discussion  of  acceptable  beam  energy  spreads  in  FELs. 
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rV.  Deviations  in  the  Relative  Phase 

To  quantify  how  the  parallel  energy  variation  affects  FEL  gain,  it  is  necessary  to 
consider  the  relative  phase  V’  of  the  electrons  in  the  ponderomotive  wave, 

dV’/dz  =  k  +  —  u)/{c/3z).  (7) 

As  is  discussed  below,  the  FEL  gain  is  directly  determined  by  the  behavior  of  the  relative 
phase  i).  In  the  small  signal  limit  (ur  — ►  0,  where  an  is  the  normalized  radiation  field),  the 
electron  energy  is  constant  and  the  axial  electron  velocity  is  given  by  /?*  =  /3q  where  0o 
is  the  initial  normalized  electron  velocity.  It  is  convenient  to  write  where 

/3xo  is  the  ideal  electron  wiggle  velocity  in  the  absence  of  field  errors  and  \80±l  0±_o\^  <<  1. 
The  deviation  in  phase  Stj)  due  to  the  field  errors®"®  is  given  by 

6^  =  -y/Io  ■  80,,  +  80l),  (8) 

where  =  0l  —  The  specific  behavior  of  Stj}  depends  on  the  specific  behavior  of 

the  transverse  motion  80x  arising  from  the  field  errors.  This  motion  has  been  examined 
in  detail  in  Ref.  6  and  is  summarized  in  Appendix  A.  In  the  following,  the  mean  phase 
deviation  {6ij;)  will  be  determined  for  (i)  helical  wigglers,  (ii)  linear  wigglers  with  flat  pole 
faces,  (iii)  linear  wigglers  with  parabolic  pole  faces  and  (iv)  an  average  of  {8xp)  over  a 
wiggler  period  will  also  be  determined. 

(i)  Helical  wigglers.  Consider  a  helical  wiggler  (with  weak  focusing)  described  by 
the  normalized  vector  potentials  given  by  Eq.  (Al)  in  Appendix  A.  The  deviation  in  the 
transverse  electron  motion  arising  from  the  error  Sa±  is  given  by 

80Z  =  {a„/2'y)kl,6y^  cos  +  60^ ,  (9a) 

60y  =  /2'f)kl,8x^  sin  k^z  +  60y,  (96) 

where  the  orbit  deviations  8x,  8y,  60^  and  80y  are  given  by  Eq.  (A3)-(A6).  Statistically 
averaging  the  phase  deviation  818  over  an  ensemble  of  wigglers  gives 


6 


+  ((^-S*)  -  sin2fc„2  +  ^  (cos2fc^2  -  1) 


X. 

^  ‘i(fcj,  -  k%)  + 


1  /  cos2(fc^  +  k0)z  cos2(fc^  —  fc^)z 


Skff  ( 


(ASjo  + 


(^to 


)]} 


(ii)  Flat  pole  faces.  Consider  a  linear  wiggler  with  flat  pole  faces  described  by  the 
normalized  vector  potential  given  by  Eq.  (All).  The  deviation  in  the  transverse  electron 
motion  arising  from  the  field  errors  is  given  by 


=  (a„/2'y)fc^iy^cosfc.„2  +60.. 
60.  =  60., 


(11a) 

(lit) 


where  the  orbit  deviations  6y,  6^y  and  60^  are  given  by  Eqs.  (A4),  (A6)  and  (A13). 
Statistically  averaging  the  phase  deviation  Srf;  gives 


r  2  1 

+  {6Bl)  sm2ky,z  + {cos2k^z  ~  1) 


^  ^(ki-kl)^  Skff  V  {k^  +  kp) 


1  /  cos  2{ky,  +  k0)z  cos  2{k^  —  k0)z 

'^0  \  (^tU  ^0)  i^w  ~  ^0) 


)]}• 


(iiij  Parabolic  pole  faces.  Consider  a  linear  wiggler  with  parabolic  pole  faces  described 
by  the  normedized  vector  potential  given  by  Eq.  (A16).  The  deviation  in  the  transverse 
electron  motion  arising  from  the  field  errors  is  given  by 


601  -  (a^/47)fc^(6x^  +  6y^)cosky^z  +  60^, 
60y  =  -{a^/2'y)kl,6x6y  cosky,  z  +  60y, 


(13a) 

(136) 


where  the  orbit  deviations  6x,  6y,  60^  and  b0y  are  given  by  Eq.  (A3)-(A6)  with  k^  replaced 
by  kff/  \/2.  Statistically  averaging  the  pheise  deviation  over  an  ensemble  of  wigglers  gives 


2ky,  Akl  1)  +  2(2ifc2  _  fc2) 

1  cos(2fc„  +  \/2k0)z  cos(2t^  —  y/2k0)z  1 

{y/iky,  +  k0)  {V^ky,  -  fc^)  .  J 
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(iv)  Wiggler  averaged  result.  Notice  that  if  the  above  results  for  the  mean  phase 
deviation,  Eqs.  (10),  (12)  and  (14),  are  averaged  over  a  wiggler  period,  then  to  leading 
order  (^V’)  is  given  by 

w)  =  {(‘bd  +  (««;>)  (15) 

where  the  resonance  condition  u/Jc  =  has  been  used  and  =  Zcy  has 

been  assumed.  Here,  1  ^  helical  wiggler  and  7^^  =  1  +  a^/2  for  a  linear 

wiggler.  Equation  (15)  is  simply  the  result  for  the  phase  deviation  as  obtained  from  ID 
theory  in  which  transverse  gradients  (weak  focusing)  are  neglected,  i.e.  kp  =  0.  Hence, 
it  is  clear  that  transverse  weak  focusing  (finite  fc^)  does  not  significantly  reduce  the  mean 
phase  deviation.  (It  should  be  mentioned  tha  in  the  trapped  particle  regime,  the  effects 
of  the  synchrotron  motion  of  the  electrons  may  reduce^ 

Physically,  niay  be  interpreted  as  an  oscillation  of  the  ponderomotive  well  due  to 
field  errors.  Maintaining  FEL  resonance  requires  to  be  small  compared  to  27r,  i.e.,  the 
width  of  the  well.  In  the  low-gain  regime,  this  phase  deviation  must  be  kept  small  over  the 
entire  wiggler  length  L.  Requiring  |  {8tl}{z  =  L))  |  <<  27r  implies  SBrm$  <  !/(’'■  •^)  ~  0.3% 
for  N  =  100  (where  ~  1  has  been  assumed  and  SBrmt  =  =  (SBy)^^^).  This  is 

the  same  condition  as  obtained  above  from  considering  the  effective  energy  spread.  In  the 
high-gain  regime,^  the  situation  is  somewhat  different,  since  the  length  scale  over  which 
the  FEL  resonant  interaction  occurs  is  the  e-folding  length  l/F,  where  F  is  the  spatial 
growth  rate  of  the  radiation.  Maintaining  resonance  in  the  high-gain  regime  corresponds 
to  keeping  small  over  an  e-folding  length:  [  =  l/E))  1  <  Since,  typically 

1/r  <<  L,  one  expects  the  high-gain  not  to  be  strongly  affected^  by  the  phase  deviation 
(in  contrast  to  the  low-gain). 
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V.  Degradation  of  FEL  Gain 

In  principle,  the  magnitude  of  the  random  walk  of  the  electron  beam  centroid  may 
be  kept  as  small  as  desired  through  the  use  of  transverse  focusing.  This,  however,  is  not 
the  case  with  the  phase  deviation,  as  is  discussed  in  the  previous  section.  In  the  following, 
the  effect  of  the  phase  deviation  on  the  FEL  gain  in  the  low-gain  regime  is  examined 
quantitatively.  In  determining  the  FEL  gain,  a  number  of  assumptions  are  made.  It  is 
assumed  that  overlap  is  maintained  between  the  radiation  and  electron  beam,  i.e.,  the 
ran'-Sm  walk  of  the  electron  beam  centroid  remains  smaller  than  the  beam  radius.  Also, 
since  weak  focusing  (transverse  gradients)  is  ineffective  in  reducing  the  phase  deviation, 
the  gain  wi'l  be  considered  in  the  ID  limit.  The  effects  of  tapering  are  neglected  and  a 
non  -diffracting,  plane  wave  radiation  field  is  assumed.  Furthermore,  the  effects  of  coupling 
to  higher  order  harmonics  (for  linear  wigglers)  will  be  neglected.  For  a  relativistic  electron 
beam,  the  normalized  amplitude  gain,  G,  is  related  to  the  relative  phase  of  the  electrons, 
V',  by’^ 

G  =  f  dz'fsin^lo,  (16) 

Jo 

where  [...)o  signifies  an  average  over  the  intial  phase  of  the  electrons.  The  quantity  G  is 
proportional  to  the  standard  definition  of  the  small  signal  gain.^^ 

The  relative  phase  V’  may  be  determined  in  the  small-signal  regime  for  which  >> 
6a\  »  aj,.  The  relative  phase  may  be  written  as  ^  +  rl>\ ,  where  is  the  relative 

phase  in  the  absence  of  the  radiation  field, 

=  V’o  +  -I-  ^^’(z),  (17) 

and  where  V'l  is  the  phase  contribution  resulting  from  the  radiation  field.  Here,  V’o  is  the 
initial  phase  of  the  electrons,  /i  is  the  normalized  frequency  mismatch, 

/i  =  -(w  -  u;o)/wo,  (18) 

where  u/q  =  cJb„(l  +  0io)Pzol^ Is  the  resonant  frequency,  and  6y()  is  the  phase  deviation 
due  to  random  field  errors  as  is  pven  by  Eq.  (8).  The  phase  contribution  resulting  from 
the  radiation  field,  is  determined  from  the  pendulum  equation^^ 

=  ^a^afisinV’^®^(2),  (19) 

dz^  71 
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which  gives 


(20) 


AL.2  /■* 

^1(2)  =  I  dz’(z  —  z')sinil}^°\z'). 

7j.  Jo 

Assuming  |  <<  1,  the  expression  for  G  may  be  expanded  giving 

G  =  l  d^’{^  1  cos  (21) 

Inserting  the  expression  for  il)i,  Eq.  (20),  indicates  that  *he  normalized  amplitude  gain  is 
proportional  to  G,  where 

G  =  dz'  J‘  dz"  [(2'  -  2") sin  (V»^®)(2')  -  ’  (22) 

where  [sin(V>^**H^')  +  ^  has  been  used.  Averaging  over  an  ensemble  of  wig- 

glers  gives  an  expression  for  the  mean  normalized  gain,  (G), 

(G)  =  f  dz'  f  dz"  {z' —  z")  {sin  [fikw{z' —  z")  +  ASrI)  (z' ,z")]) ,  (23) 

Jo  Jo 

where  ASilf(z',z"}  =  Si/>(z')  —  Si}{z").  When  =  0,  (G)  gives  the  normalized  mean 
gain  in  the  absence  of  field  errors. 

The  precise  evaluation  of  the  ensemble  average  (G)  is  nontrivial.  In  particular,  the 
evaluation  of  (G)  is  dependent  upon  the  statistical  distribution  of  the  phase  deviation  Sil>. 
Consider  random  field  errors  Sax  which  are  Gaussian  distributed.  If  the  terms  linear  in  the 
field  errors  (terms  proportional  to  SPx)  dominate  the  integrand  in  the  expression  for  Stf;, 
Eq.  (8),  then  will  tend  to  be  Gaussian  distributed.  If  the  terms  quadratic  in  the  field 
errors  (terms  proportional  to  S0\)  dominate  the  integrand  in  the  expression  for  Sij},  Eq. 
(8),  then  Stj}  will  tend  to  be  Gamma  distributed.  More  specifically,  the  statistical  behavior 
of  St/}  depends  on  whether  the  linear  or  quadratic  terms  dominate  in  the  expression  for 
variance  a  of  S^,  where  tr^  =  {S-^^)  -  {S^})^ .  It  may  be  shown  (see  Appendix  B)  that 
for  long  wigglers,  in  which  the  effects  of  field  errors  are  important,  <t  is  dominated  by  the 
quadratic  terms  in  Eq.  (8).  Hence,  the  linear  terms  in  Eq.  (8)  may  be  neglected  and  the 
phase  deviation,  in  ID,  may  be  approximated  by 

= -(*:«/7i)  r  dz'Sa\{z').  (24) 

Jo 
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Furthermore,  the  variable  is  assumed  to  be  Gamma  distributed.  Hence,  the  statistical 
average  {G)  may  be  evaluated  using  the  Rice-Mandel  approximation,^’**  giving 

(G)  =  + 

X  sin  [tiky,{z'  -  z")  +  /tan~*  {{ASrJ})//)]  , 

where  /  =  {ASiff)^ /  ~  A  more  detailed  discussion  of  the  statistical 

evaluation  of  (G)  is  given  in  Appendix  B.  Equation  (25)  describes  the  reduction  of  the 
mean  gain  due  to  random  fields  errors  in  the  low— gain  regime  for  an  untapered  wiggler. 

The  mean  of  the  quantity  ASi/j,  as  well  as  the  square  of  its  variance,  may  be  calculated 
analytically.  Using  Eq.  (24),  one  finds 

(A<V.)  =  ((«fl=)  +  (iBj))  r,  -  z"’)  ,  (26) 


{A«V’)  -  (A<V.)’  =  ^  +  (««;))%’  (z'  -  z")’ 
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X  (z'’+3z''’  +  2zV'), 


(27) 


where  Ze  =  Zex  =  Zey 


The  mean  gain  {G)  (normalized  to  the  maximum  gain  in  the  absence  of  field  errors) 
is  a  function  of  only  two  parameters:  the  product  of  the  frequency  mismatch  with  the 
number  of  periods,  fiN,  and  the  mean  phase  deviation  at  the  wiggler  end. 


=  (S’I’U  =  i))  =  -  ^  {{fBl)  +  {6Bl))  z,L\ 

Using  these  two  parameters,  (G)  may  be  written  as 

(G)  =  L'  f'dz'  f‘  dz"{z'  -  z”)(l  + 

Jo  Jo 

X  sin  [2Tr fiN{z'  ~  z”)  +  /tan“*  ((^V')ma*fi)]  , 

where  i  =  z/L  and 

/  =  (3/4)  (i'  +  z"f  /  (i'2  +  31"^  +  i'i")  , 
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(28) 


(29) 


(30a) 

(305) 


Equation  (29)  indicates  that  (G)  decreases  as  {8ilf)max  increases.  In  a  similar  fashion,  it 
is  possible  to  calculate  an  expression  for  the  variance  of  the  gain,  the  result  of  which  is 
given  in  Appendix  D.  This  variance  tends  to  be  large,  as  is  indicated  by  the  numerical 
simulations  discussed  below. 

Equation  (25)  may  be  evaluated  numerically  to  determine  the  behavior  of  the  mean 
gain.  Figure  1  illustrates  this  behavior,  in  which  the  mean  gain  [G)  is  plotted  as  a  function 
of  the  frequency  mismatch  parameter  fiN  for  several  values  of  normalized  rms  field  error 
8Brms  (0.0%,  0.1%,. ..,0.5%).  The  parameters  in  Fig.  1  correspond  to  a  linearly  polarized 
wiggler  with  B.u,  =  5.4  kG,  =  2.8  cm,  N  =  130  and  7  =  350  in  the  limit  kp  =  Q 
(transverse  focusing  is  neglected).  Notice  that  with  increasing  rms  field  error,  the  maximum 
mean  gain  decreases  and  the  position  of  this  maximum  moves  to  higher  values  of  frequency 
mismatch.  Figure  2  shows  the  peak  gain  {G)m<ix  ^  a  function  of  normalized  rms  field  error 
8Brm.$i  as  obtained  from  Eq.  (25),  for  the  above  parameters.  The  x ’s  in  Fig.  2  are  the  result 
of  an  FEL  simulation  code  for  individual  wiggler  realizations  (particular  arrangements  of 
random  field  errors).  In  these  simulation  runs,  a  random  field  error  model  similar  to  that 
of  Kincaid*'*"*  was  used  along  with  an  electron  beam  of  current  2.0  A  with  an  emittance 
of  10  /im-rad.  Notice  that  the  large  spread  in  the  simulation  results  indicates  a  relatively 
large  variance  of  the  gun. 

It  is  also  possible  to  calculate  the  effect  of  wiggler  errors  on  the  spatial  growth  rate  in 
the  high-gain  regime.*  Numerical  results  (for  a  linear  wiggler  with  B^  =  2.4  kG,  A„  =  8.0 
cm  and  X  =  15  m;  and  an  electron  beam  with  energy  50  MeV,  current  1.5  kA  and  emittance 
4.4  mm-mrad)  indicate  that  even  for  large  normalized  rms  field  errors,  SBrmt  =  0.5%,  the 
mean  spatial  growth  rate  is  only  slightly  reduced  (by  <  4%).  This  is  in  agreement  with 
the  discussion  presented  at  the  end  of  Section  IV. 
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VI.  Beam  Steering 


One  method  for  reducing  the  detrimental  effects  of  field  errors  is  through  the  use 
of  beam  steering^~**’  (external  fields  are  used  to  steer  the  electron  beam  back  to  axis). 
Analytically,  this  may  be  modeled  by  injecting  the  electron  beam  with  an  initial  transverse 
velocity  such  that  the  centroid  displacement  is  zero  at  the  end  of  the  wiggler  8x{z  = 
L)  =  0.  The  intial  transverse  velocity  may  be  specified  in  terms  of  the  perturbed  transverse 
velocity  in  the  absence  of  steering  S/3±isi  by  the  relation 

=  iz'l0i.NW),  (31) 

where  6/3±n  is  given  by  Eq.  (2).  In  the  ID  limit,  6Pxn{z)  =  8a±{z)l'y. 

Using  the  above  expression  for  one  may  calculate  the  electron  motion  in  the 
presence  of  the  field  errors  including  the  effects  of  beam  steering.  For  example,  the  mean 
square  transverse  orbit  deviation  in  the  absence  of  transverse  focusing  [kp  =  0)  is  given  by 


where  c  =  1  with  steering  and  e  =  0  without  steering,  as  derived  previously  by  Kinceud.* 
Notice  that  with  steering,  the  rms  transverse  orbit  displacement  is  maximum  at  z  =  i/2 
and  is  equal  to  1  /4  the  value  of  the  rms  displacement  obtained  at  the  end  of  the  wiggler 
in  the  absence  of  steering,  i.e., 

{8x^{z  =  1/2,  e  =  =  (l/4){6x^(z  =  L,c  =  (33) 

Similarly,  the  phase  deviation  with  (c  =  1)  and  without  (c  =  0)  steering  in  the  absence 
of  transverse  focusing  {kp  =  0)  is  given  by 

(«V,)  =  -^  z.  |r’  +  (l  -  ^  +  ^)J  .  (34) 

where  it  has  been  assumed  Zc  =  Zcc  =  Both  with  and  without  steering,  the  mean 
phase  deviation  reaches  a  maximum  at  z  =  i.  In  particular,  notice  that  the  effect  of 
steering  is  to  reduce  the  mean  phue  deviation  at  the  wiggler  end  by  a  factor  of  1/3, 

{S^{z  =  L,c  =  1))  =  (1/3)(M^  =  0)).  (35) 
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It  is  also  possible  to  calculate  the  effect  of  steering  on  the  mean  gain,  {G).  Again, 
Eq.  (25)  applies,  where  the  mean  and  variance  of  the  quantity  A^V’  including  the  effects 
of  steering  are  given  in  Appendix  C.  Similarly,  an  expression  for  the  variance  of  the  gain 
including  steering  is  given  in  Appendix  D. 

The  effect  of  beam  steering  at  the  wiggler  entrance  on  the  phase  deviation  Sip  is 
illustrated  in  Fig.  3  for  the  cases  (a)  without  steering  and  (b)  with  steering.  Here  the  solid 
curves  represent  the  mean  {Sip)  and  the  dashed  curves  represent  one  standard  deviation 
about  the  mean  (Sip)  ±  or,  where  er  is  the  variance  of  the  phase  deviation.  These  plots 
are  for  a  linearly  polarized  wiggler  with  =  5.4  kG,  =  2.8  cm,  N  =  130,  7  =  350 
and  SErms  =  0.3%  in  the  limit  kp  =  0  (transverse  focusing  is  neglected).  Notice  that  the 
effect  of  steering  at  the  wiggler  entrance  reduces  {Sip)  by  1/3  at  the  end  of  the  wiggler, 
as  is  indicated  by  Eq.  (35).  Also,  notice  that  steering  has  reduced  the  variance  of  the 
phase  deviation  by  an  equally  significant  amount.  For  cases  in  which  kp  ^  0,  it  is  possible 
to  show*  that  steering  reduces  the  mean  phase  deviation  when  the  length  over  which  the 
steering  in  performed  is  less  than  the  betatron  wavelength,  L,  <  \p.  For  cases  in  which 
Lt  >  ^0,  beam  steering  may  increase  the  value  of  {Sip). 

The  effect  of  beam  steering  at  the  wiggler  entrance  on  the  FEL  gain  (in  the  low- 
gain  regime)  is  illustrated  in  Figs.  4-6.  In  Fig.  4,  the  mean  gain  {G)  including  the  effects 
of  steering  is  plotted  as  a  function  of  the  frequency  mismatch  parameter  fiN  for  several 
values  of  normalized  rms  field  error  SErmi  (0.0%,  0.1%,.. .,0.5%).  The  parameters  in  Fig. 
4  correspond  to  a  linearly  polarized  wiggler  with  Byj,  =  5.4  kG,  =  2.8  cm,  N  =  130  and 
7  =  350  in  the  limit  kp  =  0  (transverse  focusing  is  neglected).  Figure  4,  in  comparison  to 
Fig.  1,  clearly  indicates  that  the  mean  gain  is  enhanced  through  the  use  of  steering.  For 
example,  for  SEmt  =  0.3%,  steering  increases  the  peak  gain  by  a  factor  of  approximately 
2.5.  Figure  5  illustrates  this  comparison,  in  which  the  peak  gain  {G)maxi  with  and  without 
the  effects  of  steering,  is  plotted  as  a  function  of  normalized  rms  field  error  SErmt  for  the 
above  parameters.  In  Fig.  6,  the  peak  normalized  gain  {G)max  is  plotted  as  a  function  of 
the  maximum  mean  phase  deviation,  {Sip) maxi  for  the  above  parameters  without  steering 
in  Fig.  6(a)  and  with  steering  in  Fig.  6(b).  The  curves  in  Fig.  6  remain  unchanged  for 
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various  values  of  SBrmt  and  hence,  the  maximum  normalized  mean  gain  (G)mor  is  a 
function  of  only  To  avoid  significant  reductions  in  the  mean  gain,  Fig.  6  indicates 

that  it  is  necessary  to  have  \{6il))\  «  2ir.  Figure  7  shows  the  peak  mean  gdn  {G)maz, 
including  the  effects  of  steering,  plotted  as  a  function  of  the  normalized  rms  field  error 
for  the  above  parameters.  Included  in  Fig.  7  is  the  variance  of  the  normalized  gain, 
as  obtained  from  Appendix  D,  for  several  values  of  SEmt-  A.s  was  indicated  by  the  kinetic 
simulations  for  individual  wiggler  realizations,  the  variance  of  the  gain  tends  to  be  large 
and  increases  with  increasing  rms  field  error. 


VII.  Error  Reduction  Techniques 


Several  methods  exist  for  reducing  the  detrimental  effects  of  wiggler  errors.  Above 
it  was  discussed  how  steering*"^®  the  electron  beam  at  the  entrance  of  the  wiggler  may 
improve  FEL  performance.  This  concept  may  be  generalized  to  the  case  of  multiple  beam 
steering, in  which  the  electron  beam  is  steered  back  to  axis  in  several  places  along 
the  length  of  the  wiggler.  It  may  be  shown  (see,  for  example.  Ref.  8)  that  in  order  to 
reduce  the  phase  deviation,  it  is  necessary  to  perform  steering  over  segments  of  length  Ls 
shorter  than  a  betatron  wavelength,  Ls  <  In  addition  to  beam  steering,  one  may 
consider  wiggler  errors  which  are  correlated.®  The  results  discussed  above  are  for  wigglers 
with  random  errors  which  are  assumed  to  be  uncorrelated  for  separation  distances  greater 
than  Zc  A^/2.  By  considering  a  wiggler  in  which  the  error  for  a  given  magnet  pole  is 
correlated  to  the  errors  of  the  surrounding  poles,  one  may  construct  beneficial  correlations 
which  reduce  the  detrimental  effects  of  the  errors.® 

Alternatively,  one  may  reduce  the  detrimental  effects  of  the  errors  by  considering 
an  optimal  arrangement  of  the  magnet  poles. That  is,  the  magnet  poles  are  to  be 
arranged  in  such  a  way  that  the  detrimental  effects  of  the  error  of  a  given  pole  tend  to 
cancel  those  of  the  surrounding  poles.  More  specifically,  the  magnet  poles  are  arranged 
in  such  a  way  as  to  minimize  an  appropriate  “cost  function”.  For  example,  one  may 
choose  to  arrange  the  poles  such  that  the  magnitude  of  random  walk  |^x|  is  minimized, 
where  ~  /  dz'  —  z)SBy{z').  (Notice  that  minimization  of  \  J  dzSB\  does  not 

correspond  to  minimization  of  |^®|.)  However,  the  results  discussed  above  indicate  that  a 
more  appropriate  cost  function  for  low-gain  FELs  is  the  magnitude  of  the  phase  deviation 
~  /  dz'{20xfi^$s.  +  as  is  given  by  Eq.  (8).  By  minimizing  one  reduces 

the  amount  of  gain  loss.  Ideally,  one  would  like  to  maximize  the  actual  expression  for  the 
gain,  Eq.  (22),  but  the  functional  dependence  of  the  gain  on  the  field  errors  appears  much 
too  complicated  to  be  of  practical  usefulness. 
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VIII.  Conclusions 


An  electron  beam  traveling  through  a  magnetic  wiggler  with  finite  field  errors  experi¬ 
ences  random  v,  x  6Bx  forces  which  perturbs  the  beam  motion.  This  leads  to  a  random 
walk  of  the  beam  centroid,  Sx,  as  well  as  a  random  deviation  in  the  relative  phase  of 
the  electrons  in  the  ponderomotive  wave,  6i{).  In  principle,  the  transverse  displacement  of 
the  beam  centroid,  Sx,  may  be  kept  as  small  as  desired  through  the  combined  effects  of 
transverse  beam  focusing  and  beam  steering.  This,  however,  is  not  the  case  for  the  phase 
deviation  Sij}.  Transverse  beam  focusing  is  found  to  be  ineffective  in  reducing  the  mean 
phase  deviation  {{Sijf)  is  independent  of  fc^).  Beam  steering*  may  be  used  to  reduce  |^V’I 
only  when  L5  <  A^.  As  an  example,  for  the  case  kp  =  0  and  using  steering  at  the  wiggler 
entrance  indicates  that  the  mean  phase  deviation  at  the  wiggler  end  is  reduced  by  a  factor 
of  1/3.  The  phase  deviation  leads  to  a  reduction  of  PEL  gain  (the  low-geiin  regime  is 
affected  more  strongly  than  the  high-gain  regime).  The  normalized  mean  g^un  was  calcu¬ 
lated  and  found  to  be  a  function  of  only  two  parameters,  fiN  and  (^V')mou  as  indicated 
by  Eq.  (29).  To  avoid  significant  loss  of  gain  in  the  low-gain  regime,  it  is  desirable  to  keep 
\Sil>\  <<  27r.  In  particular,  requiring  |(^V’)1  <<  2ir  gives,  using  Eq.  (15), 

{iSBD  +  iSBl))''"  (36) 

where  a  =  (1  -I-  for  a  helical  wiggler  and  a  =  (1  -1-  a^/2)^/*/o„  for  a  planar 

wiggler.  For  example,  a  helical  wiggler  with  N  =  100,  —  1  and  SBrmt  =  {SB^y^^  = 

{SByY/^  implies  that  the  normalized  rms  field  error  must  satisfy  SBrmi  <  0.3%.  Possible 
error  reduction  techniques  include  multiple  beam  steering,*’^*  correlation  of  field  errors* 
and  optimal  arrangement  of  magnet  poles. To  reduce  the  detrimental  effects  of  field 
errors  on  the  gain  in  low-gain  FELs,  an  optimal  arrangement  of  poles  corresponds  to 
minimization  of  16^1,  where  is  given  by  Eq.  (8). 
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Appendix  A: 

Transverse  Orbit  Deviations 


This  appendix  summarizes  the  results  of  Ref.  6  in  which  the  transverse  orbit  deviations 
etrising  from  random  wiggler  field  errors  were  calculated  for  (i)  helical  wigglers,  (ii)  linear 
wigglers  with  flat  pole  faces  and  (iii)  linear  wigglers  with  parabolic  pole  faces. 

(i)  Helical  Wigglers.  Consider  a  helical  wiggler  (with  weak  focusing)  described  by  the 
normalized  vector  potential 

o*  =  a«,(l  +  kl,y^  /2)cosky,z  +  ia*(z), 

(Al) 

a„  =  +  kl,x  / 2)  sin  ky,z  ->r  Say(z), 

where  it  is  assumed  k\,x^  <<  1  and  k\,y^  <<  1.  Here,  Sa^  and  Say  are  related  to  the  field 
errors  SBt  and  SBy  by 


^Uj(z)  —  a^jjk^  j  S/Z  SBy^z  ) / H,!, , 

Jo 

Say(z)  =  —ay,k^  I  dz'SBx{z')/By,. 
Jo 


The  deviations  in  the  transverse  electron  orbit  arising  from  the  errors  Sax  we  given  by 


Sx{z)  =  —  — f  dz*  sink p{z'  —  z) 

Jo 

Sy(z)  =  - ”  -  f  dz' sin  kp{z'  —  z)  — 

'ykp  Jff 


..SByjz') 
^  B^  ' 
SB,{z') 
B^  ’ 


which  correspond  to  the  normalized  transverse  velocity  deviations 


S0riz)  =  —  r  dz\oskp{z'  -  r)- 
7  Jo 

S0y{z)  =  f  dz'  cos  kp{z'  -  z) 

7  Jo 


SByjz') 
By,  ’ 
JB,{z') 

^  By,  ' 


where  kp  =  ay,ky,/y/2‘y.  Statistically  averaging  over  an  ensemble  of  wigglers  gives  the 
meam  square  quantities 


sin  2kpz' 
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{8y^)  =  1 

(a^k^y  /SBl\ 

\  1/h  J  \Bl  / 

/  sin2fcflz\ 

r  2kp  )  ’ 

(A8) 

m)  =  1 

(a^k^y  /6Bl\ 
\  7  )  \Bl, 

/  sin2li;Az\ 

r  2*;)- 

(A9) 

{sfiD  = 

/ / 6Bl^ 
\  7  )  \BIj 

V  ^c* 

'  2 

/  sin  2fcflZ  \ 

(AlO) 

(ii)  Flat  pole  faces.  Consider  a  linear  wiggler  with  flat  pole  faces  described  by  the 
normalized  vector  potential 


a  =  cosh  k^y  cos  +  6a*(z)e*  +  Say{z)ey. 


{All) 


The  X  component  of  the  orbit  deviation  is  described  by 

a„ib,„  r  .  .  r  ._„8By{z") 


8x{z)  = 


7 

a^  Aji 


/  dz'  r  dz' 
Jo  Jo 

--f 


,8By{z>) 

B^  ’ 


(^12) 

(^13) 


and  the  y  component  of  the  orbit  deviation  is  described  by  Eqs.  (A4)  and  ( A6).  Statistically 
averaging  over  an  ensemble  of  wigglers  gives  the  mean  square  quantities 


/^£iV 


m) 


7  ) 


\Bl/ 


ZcyZ, 


(414) 

(415) 


where  (6y*)  and  {8By)  are  given  by  Eqs.  (A8)  and  (AlO). 

( iii )  Parabolic  Pole  Faces,  Consider  a  linear  wiggler  with  parabolic  pole  faces  described 
by  the  normalized  vector  potential 


o*  =  co8h(fc„x/\/5)cosh(fc„y/v^)cos  k„z  +  8at{z), 

(A16) 

Gy  =  —  sinh{ky,x/y/2)  smh{k„y/ Vi)  COS  k^z  +  8ay{z). 

The  deviations  in  the  transverse  electron  orbit  are  given  by  Eqs.  (A3)-(  A6)  with  kp  replaced 
by  kp/V2.  Similarly,  the  mean  square  quantities  are  given  by  Eqs.  (A7)-(A10)  with  kp 
replaced  by  kp/V^. 
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Appendix  B: 

The  Rice-Mandel  Approximation 

Evaluation  of  the  average  (G),  given  by  Eq.  (23),  is  dependent  upon  the  statistical 
distribution  of  the  phase  deviation  Siif.  The  expression  for  6i{)  consists  of  terms  linear  in 
the  field  error  Sa±  as  well  as  terms  quadratic  in  the  field  error.  The  statistical  behavior  of 
Si/f  depends  on  whether  the  linear  or  quadratic  terms  dominate  in  the  expression  for  the 
variance  <r  of  Si/f,  where  It  is  possible  to  calculate  the  contribution 

of  the  linear  terms  to  the  variance  of  the  phase  deviation,  tri,  as  well  as  the  contribution 
of  the  quadratic  terms  to  the  variance,  ffq,  in  the  ID  limit  with  or  without  the  effects  of 
steering.  One  finds  that  the  relative  magnitude  of  quadratic  terms  to  the  linear  terms  in 
the  phase  variance  is  given  by 

4  =:  (SI) 

where  N  =  L/X^  is  the  number  of  wiggler  periods,  L  is  the  wiggler  length,  z  =  z/L  and 
2c  —  ^ui/2  was  assumed.  In  the  absence  of  steering,  g{z)  f®,  whereas  with  beam  steering 
at  the  wiggler  entrance  (one  steering  segment),  5(2)  is  given  by 

g{z)  ~  42^(1  +  3i*)-^  [(2  -  82  +  132*)  +  0.2  (-482®  +  142^)] .  (B2) 

With  steering,  Eq.  (B2)  gives  ^(1/4)  ~  0.14,  ^(1/2)  ~  0.13,  <7(3/4)  =  0.12  and  <7(1)  ~  0.20. 

As  an  example,  consider  {6B\)  ~  0.3%  and  ^(2)  ~  0.1.  Require  >>  1  indicates 

that  N  >>  26.  Hence,  for  long  wigglers  the  statistical  behavior  of  iijj  is  dominated  by  the 
quadratic  terms  and  the  linear  terms  may  be  neglected  in  Eq.  (8). 

If  the  linear  terms  are  neglected  in  Eq.  (8),  then  the  equation  for  the  phase  deviation 
reduces  to  the  generic  form 

y=  (B3) 

Ja 

where  y  represents  the  phase  deviation  and  *  represents  the  random  field  error.  If  x  is 
Gaussian  distributed  with  zero  mean,  then  y  will  tend  to  obey  a  Gamma  distribution.  The 
Rice-Mandel  approximation^*  assumes  that  the  probability  distribution  for  y,  P{y)  has 
the  general  form  of  a  Gamma  distribution.  The  parameters  occurring  in  this  general  orm 
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are  determined  by  moment  matching,  i.e.,  by  requiring  J  dyP{y)  =  1,  J  dyP{y)y  =  (y) 
and  J  dyP{y)y^  =  +  (y)^,  where  the  mean  {y)  and  the  variance  <Ty  are  assumed  to  be 

known.  One  finds 

where  a  =  {y)^/tr^.  Knowing  the  distribution  P(y)  enables  various  statistical  averages  to 
be  calculated.  For  example, 

{exp{iy))  =  J  dyP(y)exp(iy)  =  exp  tan"*  .  (B5) 

Notice  that  (y)  decreases  algebraically  as  (y)^  increases  (assuming  a  to  be  roughly  con¬ 
stant).  The  imaginary  part  of  the  above  expression  was  used  to  calculate  (G)  in  Eq. 
(25). 

As  a  final  note,  it  should  be  mentioned  that  if  the  linear  terms  in  the  expression 
for  the  phase  deviation  dominate  the  statistical  behavior,  then  y  =  J  dz'x{z')  would  be 
Gaussian  distributed.  In  particular,  (exp(iy))  =  exp(t(y)  -  <r*/2).  Hence,  if  y  is  Gaus¬ 
sian  distributed,  (exp(iy))  decreases  exponentially  as  o-*  increases.  This  implies  that  for 
Gaussian  distributed  phase  deviations  6x1),  the  mean  gain  (G)  would  decrease  much  more 
rapidly  with  increasing  than  is  predicted  by  the  Rice-Mandel  approximation,  Eq. 

(25). 
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Appendix  C: 

Statistical  Monnents  of  A6ii> 


Various  statistical  moments  of  the  function  A6ijf  =  6rl>{z')  —  8tj^{z")  may  be  calculated 
analytically  in  the  ID  limit  including  the  effects  of  beam  steering  at  wiggler  entrance 
(one  steering  segment).  In  particular,  the  mean  fbe  square  of  the  variance 

=  [ASr/?^)  —  (ASij^)^  may  be  calculated  using  Eqs.  (8)  and  (31).  One  finds 


{ASi;}  =  -  -^{SBDz^z'  -  z") 


L« 

_2  /c  d2  \2  _2  /  _/  .//\2 

<r^= - 2^(2  -2  J 


X  I  -  (z'^  +  32*'*  +  2z'z"^  ^  n  ~  +  2z' z"\ 

^  O  V  '  9  >■  ' 

-  ^  (42'  +  22"  162'*  +  92"*  +  172'^"  +  182"*2') 

9  15L  \  / 

+  ^  +  2"“)  +  19*'^"  +  242'*2"') 

where  e  =  0  without  steering  and  e  =  1  with  steering.  These  expressions  may  be  used  to 
evaluate  the  mean  gain  with  steering  using  Eq.  (25),  where  /  =  (A^V')*/<^i- 
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Appendix  D: 
Gain  Variance 


It  is  possible  to  calculate  an  approximate  expression  for  the  variance  of  the  normalized 
gain  (Toi  where 

<^1  =  (G")  -  {Gf-  {Dl) 

This  expression  is  given  here  for  completeness.  The  mean  square  normalized  gain  is  given 
by 

{G^)  =  /i  +  /2  +  /s  +  -^4  +  ht  (-^2) 

where 

h  =  f  dZi  [  dZ2  f  dz\  [  d2^A2iA22F(2i,2i,22,2i),  {D3) 

h  =  [  dzi  [  dz2  [  dz\  f  d2j^AziA22F(2i,2j,22,22).  (^4) 

Jo  J  tx  do  J  “ 

‘’"I  '‘‘'I  I  ‘'"'7,;  d2j^A2iA22^f(_)(2l,zl,22,22)>  (^5) 

74=/  d2i  f  dz2  f  dz[  f  d22^A2iA22F(_)(zi,z'l,r2,zi),  (-)6) 

Jo  Jxx  Jo  Jo  ^ 

where  A21  =  2i  —  2^  and  A22  =  22  —  zj.  The  functions  F  and  are  given  by 


F  =exp  cos(Mfc„(A2i  -  A22)  +  (A(-))] , 


2  \  -/A(i)/2 


X  cos  ftk„{Azi  ±  A22)  +  /A(±)tan  *  ^ 


(A(±))' 
/a(±)  . 


where  /a(±)  =  (A(±))Vo’i(±)  “d 


(D9) 


(A(±))  ={A(2„2;))±(A(2„2')), 


33  3 


)  +  ^Mx\.z,)  -  -  <^A(*;,x;))  • 


(7?10) 


(Dll) 
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with 


klal 


|^(*  +  y)  +  «  j  -  (* +  y)+ ^  +y* +  *y) 


(D12) 


*■«  /.< 

2  _  /  c  62  \2  ,2  /  _  ,,\2 

~X 


+  3y^  +  2xy)  +  e 


-  (5x*  -  y^  +  2®y) 


or  1 

-  ~  {Ax  +  2y-  L)  -  (l6x®  +  9y^  +  llx^y  +  18j/^®) 


(P13) 


1 


+  45^  (^*  +  24®^y^) 


}■ 


where  e  =  0  without  steering  and  c  =  1  with  steering.  Equation  (D13)  only  holds  for 
X  >  y.  When  x  <  y,  the  variables  x  and  y  need  to  be  interchanged  in  Eq.  (D13).  These 
expressions,  along  with  the  expression  for  the  mean  gain,  Eq.  (25),  may  be  used  to  evaluate 
the  gain  variance,  ffc,  with  or  without  the  effects  of  steering.  Numerically,  this  evaluation 
is  nontrivial. 
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NORMALIZED  GAIN  (G) 


FREQUENCY  MISMATCH  ixN 

Fig.  1.  Mean  gain  {G)  versus  frequency  mismatch  fiN  for  several  values  of  rms  field  error 
6Brm0  (0.0%,  0.1%, .••.0.5%)  for  a  linearly  polarized  wiggler  with  =  5.4  kG,  =  2.8 
cm,  N  =  130  and  7  =  350  in  the  limit  kp  =  0. 
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PEAK  GAIN  <G> 


5  B  rms  (%) 

Fig.  2.  Peak  mean  gain  {G)max  versus  normalized  rms  field  error  6Brm$  ior  a  linearly 
polarized  wiggler  with  =  5.4  kG,  =  2.8  cm,  N  =  130  and  7  =  350  in  the  limit 
kff  =  0.  The  solid  curve  denotes  the  theoretical  result  and  the  x ’s  denote  FEL  simulations 
for  individual  wigfder  realizations. 
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Fig.  3.  Phase  deviation  S'4>  versus  number  of  wiggler  periods  N  (a)  without  steering  and 
(b)  with  steering  for  a  linearly  polarized  wiggler  with  B^u,  =  5.4  kG,  =  2.8  cm,  7  =  350 
and  SErmi  =  0.3%  in  the  limit  kp  =  0.  The  solid  curves  represent  the  mean  {Si))  and  the 
dashed  curves  represent  one  standard  deviation  about  the  mean. 
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FREQUENCY  MISMATCH  nN 

Fig.  4.  Mean  gain  {G)  including  the  effects  of  steering  versus  frequency  mismatch  ftN  for 
several  values  of  rms  field  error  (0.0%,  0.1%,.. .,0.5%)  for  a  linearly  polarized  wiggler 

with  By,  —  5.4  kG,  A„  =  2.8  cm,  N  =  130  and  7  =  350  in  the  limit  kp  =  0. 


1.0 


NO  STEERING 


KWI, 


(b) 

STEERING 


KWLox 

Fig.  6.  Peak  gain  versus  maximum  phase  deviation  (^V’)mai  for  a  linearly  polarized 

wfggler  with  B„  =  5.4  kG,  A„  =  2.8  cm  and  7  =  350  in  the  limit  k0  =  0  for  (a)  no  steering 
and  (b)  with  steering. 


STEERING 


<5  B  rms  (%) 

Fig.  7.  Peak  gain  {G)mm*  versus  normalized  rms  field  error  SBrmi  vrith  steering  for  a 
linearly  polarized  wiggler  with  =  5.4  kG,  =  2.8  cm,  N  =  130  and  7  =  350  in  the 
limit  kff  =  0.  The  bars  denote  one  standard  deviation  as  obtained  from  the  gain  variance. 


